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We study the photo-dissociation chemistry of a quantum degenerate gas of bosonic triatomic ABC 
molecules, assuming two open rearrangement channels {AB + C or A + BC). The equations of 
motion are equivalent to those of a parametric multimode laser, resulting in an exponential buildup 
of macroscopic mode populations. By exponentially amplifying a small differential in the single- 
particle rate-coefficients, Bose stimulation leads to a nearly complete selectivity of the collective 
iV-body process, indicating a novel type of ultra-selective quantum degenerate chemistry. 



PACS numbers: 03.75.Fi, 03.75.Be, 03.75.-b 



A significant effort is currently underway with the goal 
of producing a stable Bose-Einstein condensate in a vapor 
of molecules. Several groups are pursuing schemes in 
which atomic BECs are converted to molecular BECs 
via photoassociation jjj and/or Feshbach resonance [§-§}• 
Progress is also being made towards forming a molecular 
BEC by direct evaporation of molecules in a magnetic 
trap. In this case sympathetic buffer-gas cooling is used 
to load the trap || , as opposed to the optical techniques 
used in the alkali- metal atomic vapors. 

Once the technical problems are solved, the realiza- 
tion of molecular condensates will open a new field of 
quantum-degenerate coherent chemistry. The most strik- 
ing new feature of the resulting molecular dynamics 
would be the unprecedented role of nonlinear collective 
effects, i.e. Bose enhancement and Pauli blocking, in- 
duced by the A-body quantum statistics. The proba- 
bility of any single boson to undergo a transition is en- 
hanced by the number of identical particles occupying the 
final state, whereas fermions would be subject to Pauli 
blocking of the output channels. Thus any process which 
retains a critical phase-space density would greatly dif- 
fer from that of the isolated single-particle process which 
dominates 'thermal' gas-phase chemistry j(| . 

The role of Bose stimulation in the association and dis- 
sociation of a molecular BEC has already been studied 
theoretically in the regime where a single atomic mode 
is coupled to a single molecular mode |7|-|To|. In this 
situation the entire condensate is theoretically predicted 
to undergo large amplitude coherent oscillations between 
atoms and molecules . The Bose-enhanced oscillation 
frequency is \/~NQ, where Q is the single-particle atom- 
molecule coupling frequency and N is the total number 
of condensate particles. While this prediction is signifi- 
cantly modified by the inclusion of quantum-field effects 
|TT|-|l4f , the collective evolution remains strikingly differ- 
ent from the single-pair (non-degenerate) dynamics. 

The purpose of this Letter is to extend these ideas to 
the case where more complicated molecules are dissoci- 
ated into a multiplicity of output channels. In such multi- 



channel quantum-degenerate processes there are two im- 
portant competing effects: Bose stimulation, which leads 
to exponential growth of the output channel populations, 
and condensate depletion, in which the finite number 
of initial molecules leads to competition between out- 
put channels. If certain channels have a slight advan- 
tage in coupling strength, they will deplete the conden- 
sate before other channels attain a macroscopic popula- 
tion. Related effects have been experimentally observed 
in the phenomenon of supcrradiant Rayleigh scattering 
of laser light from an elongated BEC Jl5],|l6| , indicating 
that there are no overwhelming obstacles to obtaining 
winner- takes- all scenarios in BEC systems. 




ABC 

FIG. 1. Stimulated Raman Photodissociation of ABC 
molecules 

Studying the simplest example of the dissociation of 
a triatomic ABC molecule into AB+C or A+BC pairs, 
we show that the combined effects of Bose stimulation 



1 



and competition between output channels for a finite 
initial number of molecules leads to an amplification 
of selectivity which can dramatically enhance any non- 
degenerate control scheme. Relatively small controllabil- 
ity of the single-particle process translates to near com- 
plete controllability of the Bosc-stimulated TV-body pro- 
cess. Thus, collective amplification and mode competi- 
tion effects can be employed to achieve the central goal 
of the fields of optimal- and coherent-control |l7|-|l9fl . 

We consider a BEC of triatomic ABC molecules. The 
ground-state molecules are resonantly coupled by a stim- 
ulated Raman transition to two open dissociation chan- 
nels: AB + C (AB molecules and C atoms) and A + BC 
(A atoms and BC molecules), as depicted in Fig. 1. Af- 
ter adiabatic elimination of the intermediate ABC* state 
the A-body second-quantized Hamiltonian in the interac- 
tion representation and the rotating wave approximation, 
reads 



ni^ B (r)*k(r) 



n 2 ¥ A (rW BC (r) 



h.c. 
(1) 



where J7i , ^2 are the coupling frequencies into the AB+C 
and A + BC dissociation channels, respectively and <SE>, 
$ab, ^c, ^bc are the bosonic field operators for 
ABC, AB, C, A, BC, respectively. 

Using the Hamiltonian (Q), the Heisenberg equations 
of motion for these operators are 



i— $ = n^AB^c + fy^A^BC 
at 

i~j2^ AB = ^1$% 



(2a) 
(2b) 
(2c) 
(2d) 
(2e) 



Equations (0) are closely related to the quantum equa- 



tions of motion of a multimode laser 
AB + C dissociation channel (Eqs. 

-BC channel (Eqs. 



201] , in that the 
2q)) is only cou- 
)) via the molec- 



pled to the A + BC channel (Eqs. fl2d|,N) 
ular reservoir <t. In order to demonstrate the collective 
amplification of selectivity, we use a five-mode model, ne- 
glecting the spatial variation of the various fields. This 
approximation, which assumes there is only one avail- 
able state per specie, is valid for a zero temperature 
ABC condensate at exact resonance with the dissocia- 
tion threshold of the AB + C and A + BC channels. The 
field operators of Eqs. (|^) are thus taken to represent 
the particle creation and annihilation operators of the 
respective modes. In the undepleted pump approxima- 
tion OJl4l], the intense $ field is replaced by a constant 



c-numbcr $ ~ VAf. Assuming without loss of generality 
that f2i = f2f and CI2 — ^21 the initial time evolution of 
the operators ^>ab, &A> and ^ bc is thus given as 



^ AB (t) 



^ab(0) coshwii - i& c (0) sinh^ii 
(0) coshwii + W A b(0) sinhujit 



ty A (t) = ^a(O) coshw 2 * - iW BC (0) sinhw 2 i 



W BC (t) = ^ BC (0) coshuj2t + i 1 ^ a{0) sinhu) 2 t 



(3a) 
(3b) 
(3c) 
(3d) 



where u>\ = yNQi and u>2 — y/NCl^. Assuming 
(*Ab{0)) = (*o(0)) = (*a(0)) = (*bc(0)) = (no 
initial product coherence) , the average numbers of parti- 
cles (rij) = in the modes j = AB,C,A,BC, are 
given, respectively, as 

{n>AB(t)) = (n A B{0))cosh 2 uJit 

+ (l + (ri c (0)))sinh : W , (4a) 
(no(*)> = (^c(0))cosh 2 cj 1 t 

+ (l + (nAs(0)))sinh 2 wit , (4b) 
(riA{t)) — {ua(0)) cosh 2 uo2t 

+(l + («sc(0)))sinh 2 ^ , (4c) 
(riBc(t)) = (n B c(0)) cosh 2 oj 2 t 



+(1 + (riA(0)))sinh 2 uj 2 t 



(4d) 



Starting with a pure ABC condensate (uab(Q)) — 
(nc(0)) = (ua(0)) — (n B c{0)) = 0, the average prod- 
uct populations grow as 



(5a) 
(5b) 



(n AB {t)) = (n c (t)) = sinh 2 ViV^it , 
(n A (t)) = (n BC (t)) = sinh 2 ^^ • 



The initial time evolution of Eqs. (|5|) corresponds to 
an initial nonexponential (quadratic rather than linear 
in time) spontaneous dissociation process, followed by an 
exponential stimulated process |Q . We note that exactly 
the same equations have been routinely used for over 
three decades, to describe the phenomenon of parametric 
superfluorescence pl| , p2| , i.e. the parametric amplifica- 
tion of noise photons, in quantum optics. The feature 
which produces the striking selectivity of Bose-enhanced 
chemistry is that any small difference between the two 
rate coefficients, D,2 =/= fii, is exponentially amplified by 
the stimulated dissociation, so that 



N 2 (t)/N 1 (t) « e 2VF(^-^i)*. 



(6) 



The linearized model of Eqs. (g) is only valid as long 
as the population of the triatomic mode ($t$) [ s large 
and the effect of its depletion on the products population 
growth is negligible. In order to go beyond this approx- 
imation we solve the full A-body problem numerically. 
Starting with no products and fixing the total number of 
particles N, the accessible Hilbert space is restricted to 
Fock states of the type \tiabc, tiab, nc, n A , tlbc) with 
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jiab = nc, n A = riBc and uabc = N - (n A + nc)- 
Using these states to represent the Hamiltonian (|j]) we 
solve numerically the TV-body Schrodinger equation, 



d 

dt 



n ly /N -(n 1 + n 2 ){ni + l)c ni 



+ l,n 2 



+ n$y/N-(n 1 + n 2 ) + ln 1 c* i _ li „ 2 
+ Q 2 y/N - (ni + n 2 )(n 2 + l)c nii „ 2+ i 
+ ^ y /N-{n 1 +n 2 ) + ln 1 c* niina _ 1 , (7) 

for the dynamics of the probability amplitudes 

c ni ,n 2 (t) = (TV- (m +n 2 ),ni,ni,n 2 ,n 2 |*) . (8) 
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FIG. 2. Expectation values of the AB + (?(—•—) and 

A + BC ( ) channel populations as a function of time 

for (a) N — 1 and (b) TV = 1000. The coupling rate ratio is 
fi 2 = 2Qi. 

In Fig. 2, we plot the expectation values (nc), (nA), 
corresponding to the populations in the AB + C channel 
and the A + BC channel, respectively, as a function of 
time, for Q 2 = 20 x . When TV = 1 (Fig. 2(a)) the ratio 
between the two channel populations is, as expected for 
the single-particle process (ha) / (nc) — |^W^i| 2 = 4- 
However, for TV = 1000 (Fig. 2(b)), the Bose stimulated 
population ratio is (n A )/{nc) > 40. Thus, the reaction 



outcome is highly dependent on the total number of par- 
ticles, becoming more selective as TV increases. 
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FIG. 3. Channel population ratio (tia) / (nc) when 
(nA) = ( n A) max as a function of TV for Q.2 = 1.25f2i 
(dash-dot), O2 = 1.5Sli (solid), and Q.2 = 2fii (dashes) 

Our numerical calculations are limited by available 
computation power to TV < 1000. However, the num- 
ber of particles in current BEC experiments in trapped 
dilute alkali gases is sometimes as high as 10 7 , suggesting 
a far greater selectivity enhancement. In Fig. 3 we plot 
the channel population ratio when the leading channel 
population reaches a maximum, as a function of TV for 
various values of fi 2 /^i- The curves are linear fits to the 
numerical data, demonstrating a power-law dependence 
of the selectivity as ~ TV" . The power a increases with 
the rate coefficient ratio f2 2 /Sli, as shown in Fig. 4. 

In order to analytically estimate the dependence of the 
selectivity on the total particle number, we take a per- 
turbative approach where we first obtain the time td at 
which the leading channel depletes the condensate, and 
then evaluate the channel population ratio at td- Assum- 
ing a purely exponential growth, the depletion time is 
t d = logTV/2 v / TVft 2 . Substituting into Eq. (|), we find 
that the selectivity is enhanced as TV", where 



a = 1 - Qi/fi 2 



(9) 



is 



so that log[TV 2 (trf)/TVi(trf)] w a log TV. This result 
confirmed by the numerical calculations in Fig. 3, giv- 
ing a values of 0.17, 0.32, and 0.49 for ft 2 M = 1.25, 
fi 2 /Qi = 1.5, and 2 /f2i = 2, in good agreement with 
the predicted values of 0.2, 0.33, and 0.5, respectively. 
Extrapolating to TV = 10 6 , we find population ratios that 
are greater than 13, 100, and 1200 for fi 2 /fii = 1-25, 
f2 2 /f2x = 1.5, and Q 2 /fii = 2, respectively (as compared 
to single particle ratios of 1.56, 2.25, and 4, respectively). 
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The dependence of the selectivity amplification con- 
stant a on the rate coefficient ratio f^/f^i is plotted in 
Fig. 4. The solid line depicts the estimate of Eq. (j|). 
At the limit of a high rate coefficient ratio the selectivity 
is enhanced as N, implying that all the particles would 
dissociate into the preferred channel. We note that since 
the linearized model underestimates the condensate de- 
pletion time, the numerical values of a are consistently 
lower than the theoretical estimate of Eq. (0) . 
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FIG. 4. The variation of the selectivity amplification con- 
stant q with the rate coefficient ratio Q2/Qi- The solid line 
corresponds to the estimate of Eq. (Hj), whereas the circles 
are obtained from the numerical quantum solutions 

While the five-mode model is sufficient to produce 
Bose-enhanced gain, it neglects various spatial effects. In 
particular, channel losses due to the relative translational 
motion of the products are not accounted for under the 
assumption that energy conservation forces the reaction 
products to be at rest. In fact, by imparting additional 
momentum via the laser fields, such losses could be used 
to increase selectivity by setting a threshold for ampli- 
fication. The population in each rearrangement chan- 
nel would only be stimulated if the respective pump rate 
\/Nflj is faster than the channel depletion. Thus, if the 
motion in the AB+C channel is faster than in the A+BC 
channel it will be possible to tune the laser intensity so 
that only one pathway is Bose stimulated. In particular, 
if A and B are atoms of similar mass ttia ~ Wbi an d 
C is a much lighter atom mc <C tua, tub, the A + BC 
channel will be enhanced at much lower pump intensi- 
ties than the AB + C channel. Consequently, selectivity 
could be attained even for S7i = Q2, an issue currently 
under active investigation. 

By drawing an analogy from the phenomena of BEC 
superradiance (l^Jf^l to multi-channel reactive scatter- 
ing in molecular condensates, we have shown that the 



selectivity of reactive processes carried out in a quantum- 
degenerate gas, can be greatly amplified by collective en- 
hancement effects. This selectivity is the result of inter- 
play between bose enhancement, which leads to exponen- 
tial growth, and competition between modes for a finite 
number of initial molecules. The effect was demonstrated 
using the simple model system of two-channel photodis- 
sociation, where due to a small difference between the two 
channel rate coefficients, population was shown to be di- 
rected by the TV-body stimulated process into the favored 
channel, approaching a nearly complete selectivity in the 
limit of large N. In the quantum degenerate regime, it 
is therefore only necessary to exert slight control over 
the reaction rates in order to achieve near-total control 
over the reaction output. Future work would address the 
question of complementary Pauli-blocking effects when a 
molecular BEC is dissociated into fermionic constituents. 
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